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Abstract 

We extend the analysis of the recently obtained stringy cosmological solutions induced by 
thermal and quantum effects, once space-time supcrsymmctry is spontaneously broken by 
geometrical fluxes. Cases in which more than one modulus participating in the supersym- 
metry breaking mechanism are investigated. The free energy is obtained at the full string 
level. In the intermediate cosmological region where the temperature and the supersymme- 
try breaking scale are sufficiently smaller than the Hagedorn temperature, the quantum and 
thermal corrections are under control and calculable. The reason is that the contributions 
to the effective potential of the moduli that are not participating in the supersymmetry 
breaking are exponentially suppressed. The backreaction on the initially flat background 
results in many cases into cosmological evolutions, where the dynamics of all complex struc- 
ture moduli is frozen. The solutions describe effectively a radiation dominated era, where 
thermal effects are never negligible, even if the temperature tends to zero at late times. We 
analyze several types of supersymmetry breaking patterns and examine the stability of the 
corresponding radiation era. 
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1 Introduction 



Perhaps the most natural setting for testing string theory is the cosmology of our Universe. 
By now, there is a plethora of cosmological data favoring the phenomenological model of hot 
Big Bang cosmology, where the cosmological evolution starts with a highly singular event, 
the Big Bang, followed by an initial period of rapid inflation, a very high temperature phase, 
a proportionally large amount of dark matter and dark energy [1]. Many features of this 
phenomenological model are not well understood, and the hope is that their explanation will 
arise from a fundamental theory of quantum gravity such as string theory. 

In an effort to build a concrete theoretical framework for studying cosmology, a class of 
string theory vacua, where the backreaction of both thermal and quantum effects can be 
systematically taken into account, was recently examined in [2-5]. In particular, starting 
with weakly coupled four-dimensional heterotic string models, with initial N — A or N — 2 
space-time supersymmetry [5] , and implementing the thermal and the quantum corrections 
due to the spontaneous breaking of supersymmetry, cosmological solutions are found, at 
least when the temperature T and the supersymmetry breaking scale M are sufficiently 
below the Hagedorn temperature Th- The string coupling must also be sufficiently small 
so that pcrturbativc, thermal equilibrium computations of string theory can be applied. 
In these examples, all moduli that are not involved in the breaking of supersymmetry give 
exponentially suppressed contributions to the (thermal) effective potential. This is essentially 
the underlying reason for the no-scale structure [4,5] characterizing the models. Due to this 
remarkable property, the thermal and quantum corrections are under control and calculable. 
We would like to stress here one of the merits of the one-loop string computations we perform: 
namely, the absence of both infrared and ultraviolet ambiguities leads to a well defined energy 
density and pressure [5]. 

The backreaction of the thermal and quantum corrections on the initially flat metric and on 
certain moduli fields (including the dilaton field and the geometrical moduli participating 
in the supersymmetry breaking mechanism) induces the cosmological evolution. Both the 
temperature and the supersymmetry breaking scale evolve in time inversely proportional to 
the scale factor of the Universe [5]. As the Universe expands, the system remains weakly 
coupled in the future, and so we can trust the solution at later times. However, if we extrap- 
olate the solution back in time, the temperature rises and as this approaches the Hagedorn 
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temperature, we encounter the Hagedorn instabilities of thermal string theory, indicating a 
non-trivial high temperature phase transition [6-8] . At this point the dynamics drive the sys- 
tem outside the perturbative domain, and the large backreaction of the tachyon condensates 
involved in the transition may drive the system to strong coupling. In the literature there 
are some speculative proposals concerning the nature of this early times/high temperature 
phase transition [6-12], but still an adequate quantitative description of the dynamics is 
lacking. We will assume that such a stringy phase transition occurs at very early times, and 
that the system eventually exists in a low temperature, T << Th, weakly coupled phase, 
with three large spatial directions and possibly some intermediate size spatial directions that 
generate the scale of supersymmetry breaking, M « Th., via geometrical fluxes. In this low 
temperature phase, we have good analytical control to analyze the subsequent cosmological 
evolution, using string perturbation theory techniques. 

To be more precise, we are forced to separate the cosmological evolution in four distinct 
phases, according to the value of the temperature, namely: 

(i) The very early "Big Bang" phase, where the underlying string theory degrees of freedom 
are strongly coupled. Perhaps string dualities can be applied to understand this phase and 
resolve the classical Big Bang singularity. See e.g. [10-16] and references therein. 

(m) a high temperature stringy phase, T < Th, where string oscillators and winding states 
must be properly taken into account [6-9,17]. Often, these lead to a non-geometrical struc- 
ture of the Universe, e.g. the T-fold cosmologies of [2], or even to a change of the topology 
and dimensionality of space [10, 11]. 

(iii) The third phase has features similar to that of a standard, radiation dominated Pried- 
mann cosmology. Here, the Universe has cooled down to temperatures far below Hagedorn. 
The effects of string massive states are exponentially suppressed (9(cxp(— Ms/T)). In this 
phase, the ratio of the temperature T and supersymmetry breaking scale M is fixed, both 
evolving inversely proportional to the scale factor of the Universe [4,5]. In cases with = 1 
initial supersymmetry, the behavior can be that of a cosmological constant dominated infla- 
tionary universe [3,4]. 

(iv) At lower temperatures, new phenomena such as the electroweak phase transition, QCD 
confinement and structure formation are taking place. We expect also that in this phase, 
some dynamics becoming relevant at these lower temperatures will stabilize the no-scale 
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modulus associated to the supersymmetry breaking scale [18]. 

Some interesting ideas concerning the first two stringy phases have been presented recently in 
[12] , where it was argued that the introduction of certain chemical potentials in the standard 
canonical ensemble of superstrings removes the Hagedorn instabilities. These ensembles are 

characterized by a "Temperature duality," Z(T jTu) = Z{Th/T). (See also [19] for some 
related work.) Equally interesting are the models of [11], which possess a new kind of massive 
boson-fermion degeneracy symmetry. Type 11, heterotic and orientifold models have been 
presented. Some proposals are put forward, in the framework of these theories, concerning 
the early structure of the Universe. 

Here however, we would like to examine more thoroughly the generalization of supersym- 
metry breaking in the cases where more moduli are involved, not only in heterotic string 
but also in type 11 closed string theories, working in the intermediate region [iii). It is 
interesting that this intermediate phase can be studied with high precision at the full string 
level [5], thanks to to the fact that just below the Hagedorn temperature, the theory pos- 
sesses a no-scale structure. The free energy is set by a single, overall scale, which can be 
chosen to be either the temperature or the supersymmetry breaking scale, and the rest of 
the dependence is given in terms of functions of dimensionless, complex structure like ratios. 
The backreaction on the initially fiat background induces the cosmological evolution [4,5]. 

In the string models studied in [5], only a single modulus was participating in the super- 
symmetry breaking mechanism. In this work we extend the analysis for cases where more 
geometrical moduli participate in the spontaneous breaking of supersymmetry. As we already 
remarked, all such moduli give non-exponentially suppressed contributions to the effective 
potential. We examine in more details the case where the supersymmetry breaking arises 
via geometrical fiuxes [20] induced by stringy Scherk-Schwarz [21] boundary conditions along 
two internal spatial cycles with radii denoted by R4 and R5. When the complex structure 
ratio U = -R4/-R5 is stabilized to be of order one, the supersymmetry breaking scale M will 
be proportional to the inverse of the volume modulus y/R^R^. (The whole U dependence of 
the one-loop effective potential is presented and analyzed in this work.) The results can be 
easily extended to cases involving more internal cycles. The heterotic or type II geometrical 
fluxes are introduced by utilizing the helicity and/or other internal R-symmetry charges. 
In the type II models, left-moving and right-moving R-symmetry charges can be coupled 
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symmetrically or asymmetrically to the two cycles. 

As we will see, the low energy dynamics of some models admits a solution describing a 
radiation dominated era. During the evolution, the supersymmetry breaking scale M and the 
temperature T evolve in time the same way as the inverse of the scale factor of the Universe. 
It is then important to analyze the stability for the dimensionless modulus describing the 
shape of the internal space, U = We find a rich structure of phenomena, depending 

on whether the corresponding effective potential admits a minimum, a maximum or a run- 
away behavior. Complementary results can be found in [22]. There, it is shown that these 
solutions are attr actors for the dynamics. That is, there are basins of initial conditions 
whose associated cosmological evolutions converge towards the radiation dominated era with 
stabilized complex structures. The results of [22] are important for the following reason. Due 
to the lack of quantitative control, the ambiguities of the Hagedorn transition exit can be 
parameterized in terms of initial time boundary conditions, "capping off" the cosmology at 
an early time just after the temperature has dropped below Hagedorn. It is then important 
to know that for generic initial conditions, the subsequent cosmological evolution is described 
by a solution such as the one displayed in this paper. We expect that the degrees of freedom 
describing the physics of the early times/high temperature stringy Hagedorn phase, decouple 
at later times and low temperatures. 

In this work we choose to break supersymmetry via geometrical fluxes which can be imple- 
mented consistently not only in the effective no-scale supergravity theory, but also at the 
string one-loop level. Our work is concerned with the intermediate cosmological region {iii) 
just before the electroweak phase transition. In this regime, the geometrical moduli i?„ that 
do not participate in the supersymmetry breaking mechanism (that is, moduli other than 
the dilaton and the radial moduli R4 and R5) can: 

- Either get a soft breaking mass at tree level, see for instance [21]. They are thus thermally 
fluctuating about their ground state. 

- Or, they remain massless at tree level. 

Both types of such moduli give exponentially suppressed contributions to the one-loop ef- 
fective potential when their vevs are of order unity as compared to the string scale i.e. 
Rn <^ 1/T, 1/M. In [23], it is explicitly shown that if this condition is initially satisfled 
right after the exit from the Hagedorn phase, it remains so during the whole intermediate 
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radiation dominated era (iiijQ If we turn on kinetic energies for these massless moduli, 
these give sub-dominant contributions to the effective Hubble equation, of the order 
The same is valid for Wilson line moduli. The values of Wilson lines associated to extended 
symmetry points are preferred at finite temperature since these points give rise to more 
massless degrees of freedom. It is important however to analyze and provide a mechanism 
for the stabilization of relevant moduli just after the electroweak phase transition. One of 
these is the no-scale modulus that sets the supersymmetry breaking scale and therefore the 
electroweak symmetry breaking scale via radiative corrections. In this work we do not carry 
out this later time analysis. We also expect that the use of other types of fluxes for moduli 
stabilization (which at the moment cannot be fully implemented at the string perturbative 
one-loop level) will result in similar features in the intermediate cosmological region we are 
working with, as was suggested in [4] . 

The paper is organized as follows. In Section 2, we describe the basic principles of the 
construction, clarifying the domain of validity of our analysis. We examine different classes 
heterotic and type II string vacua and implement various ways that lead to a spontaneous 
breaking of supersymmetry. We calculate the free energy at the full string level, and obtain 
the generic structure of the thermal effective potential in the intermediate region where 
T <^ Th- We use these general results to find their counterparts at zero temperature. We 
also derive in the Einstein frame the gravitational equations and the equations of motion 
for the main moduli participating in the supersymmetry breaking mechanism, including the 
dilaton field. We present solutions at zero temperature, where the time evolution of the scale 
factor and the stability of the complex structure U are analyzed. 

In Section 3, we present explicit thermal models leading to radiation dominated cosmological 
evolutions. The compatibility between the gravitational equations and the equations of mo- 
tion of the relevant moduli leads to the equation of state p = (3+n) P, where n is the number 
of internal dimensions involved in the supersymmetry breaking mechanism. In addition, we 
analyze the stability for the modulus U. We compare to the zero temperature situation, 
mainly to show the relevance of the thermal corrections to the cosmological evolution. We 
show that during the radiation era, even when the temperature is very small, thermal effects 

^To be precise, this is true as long as the initial time-derivative of i?„ does not exceed some "escape 
velocity" , above which the modulus enters in a different regime where it is spontaneously decompactified. 
fn the latter case, the system is driven into a radiation-like era in higher dimensions. 
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are never negligible. The qualitative behaviors of thermal and non-thermal evolutions are 
drastically different. This is due to the non-linear character of the gravitational and moduli 
equations. 

The final Section is devoted to conclusions and perspectives. 



We consider initially supersymmetric and weakly coupled fiat backgrounds within the context 
of four-dimensional superstrings constructed via orbifolds [24] and/or via the free fermionic 
construction [25] . By turning on certain non-trivial geometrical fiuxes, we can spontaneously 
break space-time supersymmetry [21]. The procedure that we follow involves coupling some 
of the internal lattice quantum numbers to the space-time Fermion number F and/or to 
any of the discrete R-symmetry charges Qr- This is a generalization of the Scherk-Schwarz 
mechanism to superstrings [21]. In addition, the system is put at finite temperature [6-8]. 
The temperature and the supersymmetry breaking scales arc taken to be below the Hagedorn 
temperature. Our aim is to study cosmological evolutions induced by the thermal and 
quantum corrections, as in [5], and in particular investigate how some of the geometrical 
moduli, participating in the breaking of supersymmetry, can be stabilized around a local 
minimum. 

2.1 Heterotic models 

We start with heterotic string compactifications on six-manifolds of the form S'^(i?4) x 
S^lR^) X Here, the choice of the four-manifold determines the initial amount of 

supersymmetry: A^4 = 4 for the case of T^, and A^4 = 2 for the cases of T^/Z2 orbifold and 
compactifications. Two internal cycles, associated to the X'^ and directions, have 
been singled out, since these are to be utilized to break supersymmetry spontaneously. 

We illustrate the derivation of the pressure with the simplest example. The initially super- 
symmetric string partition function is given by 



2 General setup 



z 



V3 f d^r 



4 ^(6,22) 



(2.1) 



a,b 
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where the Euchdean time direction is compactified on a circle of radius Rq. V3 is the volume 
of the three large spatial directions. The r(6,22) lattice is associated to the zero mode contri- 
bution of the internal six-manifold, along with the x or 5*0(32) right- moving lattice. 
For instance, the Eg x Eg case on gives the block 

r(6,22) = r(6,6)^$:Cf ^E^ll']'' (2.2) 

7,5 7', 5' 

where 7, 6 and 7', 6' are integers defined modulo 2. 

We wish to implement a non-zero temperature in the model. This is done by coupling the 
momentum and winding quantum numbers associated to the Euclidean time circle to the 
space-time fermion number F [6-8] [5]. The contributions of the right-moving world-sheet 
degrees of freedom to F are always even. Thus, at the level of the one-loop string partition 
function, the operation amounts to replacing the r(i^i)(i?o) lattice with 

r(i,i)(i?o) - r(i,i)e](^o) i-ro+tHo+soko ^ (2.3) 

where is a Z2-shifted lattice [26] given by [6-8] [5] 

Ti lilm = A y e-"#l(^'^'+^)+(''^+')^l' . (2.4) 
- 

* m,n 

The spontaneous breaking of space-time supersymmetry is done by coupling the two 1) 
lattices associated with the internal circles of radii R4 and to generic R-symmetry charges 
[5-8,21]. In the case of models with N4 = 4 initial supersymmetry, all such charges associated 
to the left-moving world-sheet degrees of freedom are equivalent by symmetry. Different 
choices exist involving right-moving gauge R-charges [5]. For example, consider the Eg x Eg 
models and decompose the Eg representations in terms of 50(16) ones. One can choose 
R-charges which are odd for the 5*0(16) spinorial representations and even for the other^ 

We will present a class of cases, where starting with A^4 = 4, i^g x Eg vacua, space-time 
supersymmetry is broken if we couple the X* lattices, i = 4,5, to F + Qi, where the right- 
moving charges Qi are odd for the 5*0(16) spinorial representations associated with one or 
both the Eg factors. So we replace the X'^ and lattices as follows: 

r(i,i) (Ri) ^(I'l) [fc] (^*) ^^)i^+Q^)9^+{b+U)h.+erg.h, _ (2.5) 

^In the N2 — 2 orbifold models, one can choose R-charges associated to the twisted planes to which 
the left-moving world-sheet degrees of freedom contribute as well [5] . 
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Qi, i — 4, 5, can be identified to be either 7, 7' or 7 + 7'. [7 and 7' are odd for tlie 
corresponding SO {16) spinoriaf representations.] Lj is equal to S, S' or S + S' respectively, 
as dictated by modular invariance, and ei — 0,1 depending on the modular transformation 
T — > T + 1. Under this, 

a + Qi ^ a + Qi , b + ^ a + b + Qi + Li + €i . (2.6) 

For instance, for {Qi,Li) = (7 + 7', 5 + 5') one has = 1, while for {Qi,Li) = (7, 5) 
or (7', 5'), ei vanishes. With these modifications taken into account, the one-loop string 
partition function is given by 



7 = ^3 / ^ ^ ^(4,4) 1 Sr(_\a+b+abQ ra 



a,6 ^2.7) 



7,5 7',i5' 1=0,4,5 

In this equation, we have Qq = Lq = 0. Redefining a — d + hi and b = b + ^j, and 
using the Jacobi identity [5], one obtains 

^ I iV r(4^4) [i+g.fc»'^ 

(27 



'27r)3 i^2r^/^ ^12^24 ^ [i+EiS., 



2 



X zEE^l^]'^!^']' 11^(1.1)6;] (-)^>^'^'+^^^^+^''^^ 



Ei /^i+Ei 9i+(Ei '»i)(Ei 9i) 

(2.8) 



7,5 7',i5' i 

In the large radii regime Ri » Rh, where Rh is the Hagedorn radius, the system is free of 
tachyons. The odd winding sectors, hi — 1, are exponentially suppressed. In this regime 
only the sectors hi — 0, i — 0,4,5 and go + §4 + 95 — 1 modulo 2 contribute significantly (the 
latter condition due to the fact that 9[\] vanishes identically). Furthermore, if the internal 
lattice r(4 4) moduli are kept to be of order unity, we can express the leading contributions 
as the following integral [5]: 

Si s (2.9) 



Jo ^2 ™. 



In the first line the sum over s runs over the 2^ x 504 massless boson/fermion pairs of 
the initially supersymmetric model. The contributions of massive states are exponentially 
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suppressed, of order e The integral gives the pressure in the string frame: 

'1 (2.10) 



rtii 



[(2mo + ^o)'i?g + (2m4 + ^4)'i?l + (2m5 + ~g^fR^ ' 



(2.11) 



We parameterize the various moduh as follows: 

T := ^-^= , M := / , ReTi := R4R5 = T , 

ReC/i := — = W , w := — , 

where 5" is the 4d dilaton-axion modulus, ReS" = e"^*^^ . The two supersymmetry breaking 
scales, in the Einstein frame, are the temperature T and the scale M, which for typical 
models determines the gravitino mass scale. The pressure in this frame is related to the 
string frame pressure by: 

^ Pstring = T^P{U,U) , (2.12) 



(ReS) 
with 

p{u,U) = niooPioo(M,W) + noioPoioiu,U) + nooiPooi(M, W) + niiiPiu{u,U) . (2.13) 
The coefficients Ug^g^g^ are given in terms of the supersymmetry breaking R-charges, 
nioo = 2^x 504, n,,, ^Y.^-f'^'^ , n,o, = Y.^-f'^'^ . rzni = E(-)^'^^^^'''^^^' 

s s s 

(2.14) 

while the dependence on the complex structure moduli u and U involves the shifted Eisenstein 
functions: 

2 

P-gmi'^M)^:^ XI T T3- (2-15) 

mo,m4,m5 (27flo + QoYu^ + (2m4 + QaYU"'^ + (2m5 + ^5)^^ 

2.2 Type II models 

We construct Type II models with similar thermal and supersymmetry breaking properties. 

In these examples, the internal manifold involves cither a factor for A^4 = 8 initial super- 
symmetry or a T"^ I^Li factor for N/^ = 4. Orientifolds of these lead to models with A'4 = 4 
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and N4 — 2 initial supersymmetry respectively, and include open string matter sectors. At 
weak coupling in four dimensions, these are dual to heterotic models [27] , some of which we 
considered in the previous section. Models with N4 — 2 initial supersymmetry can also be 
constructed if we start with a T^/ (Z2 x Z2) orbifold [28] . We illustrate the derivation of 
the pressure in the intermediate cosmological region, with T <^Th, for the type II A^4 = 4 
models, but the results can be generalized to the other cases. 



The A^4 = 4 partition function is 
r rfV 1 

z ^ 



H,G 



a,b a,b 

The rVZ2 part is given by [26]: 



(2.16) 



%,4)[f] = |J|yi, when iH,G) = (0,0) 

%,4)[g] = .niXl^i^ > ^^^^ iH,G)^iO,0) 



(2.17) 



The characters H, G are integers defined modulo 2. 

As usual, the finite temperature is implemented by inserting the thermal co-cycle and re- 
placing the Euclidean time lattice as follows [6-8] : 

r(i,i)(^o) - Yl ^(1,1)601(^0) (-)("+^)5o+(;>+5).o . (2.18) 
ho, go 

In contrast to the heterotic case, the contributions to the space-time fermion number F from 
both the left-moving and right-moving sectors can be odd or even. In the sequel, we denote 
by Fl the contribution of the world-sheet left-movers to the space-time fermion number and 
similarly for Fr. 

There are several ways to break the initial A^4 = 4 supersymmetry spontaneously, either 
by symmetric or asymmetric geometrical fluxes [20,21,24]. The two left- moving space-time 
supersymmetries can be broken if we couple either or both the and lattice charges to 
Fl, or to left-moving R-charges associated with the twisted planes: Fl + Qi. Also, the two 
right-moving space-time supersymmetries are broken by coupling the lattice charges to 
or to Fr + Qi. Each lattice is replaced as follows [12]: 

1)(-Ri) — > r(ii)[~*](it!j) (^ — '^[("■+Qi)9i+{b+Li)hi+gihi]ei ^_'^[{a+Qi)gi+{b+Li)hi+gihi]ei ^ (2.19) 
hi,gi 
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where Qi, Qi can be set to zero or identified with the twist charge H. Correspondingly 
Li, Lj can be set to zero or identified with the character G. Also, we have introduced the 
parameters ej, taking the values or 1, to indicate whether we couple the circle i to the 
left- or right-movers. 

In particular, we will examine 3 distinct cases where A'4 = 4 is spontaneously broken to 
A^4 = (and then thermalized) : 

-Case 1 : Two asymmetric breakings, e.g. (£4, £4) = (1,0), (£5,65) = (0,1). 

-Case 2 : One symmetric and one asymmetric breaking, e.g. (€4,64) = (1,1), (£5,65) = 
(0,1). 

-Case 3 : Two symmetric breakings, (64,64) = (1, 1), (65,65) = (1, 1). 

In addition, we will consider a case where the N4 = A supersymmetry is partially broken to 
N4 = 2. The remaining supersymmetries are then broken by thermal effects: 

-Case 1' : Two left-moving asymmetric breakings, e.g. (64,64) = (1,0), (65,65) = (1,0). 

The partition function can be written as follows: 

1/3 f d^r 1 1 



(2< 



[ d^T 1 1 X ^ jj 



H,G 

a,b a,b 

hi,gi rhi,ni 

Here (60, eo) = (1,1) and {Qo,Qo) = {LoM) = (0,0). 

As in the heterotic case, we are interested in the regime where the radii Ri, i — 0,4,5, 
are much bigger than the Hagedorn radius, Ri ^ Rh- In this intermediate cosmological 
regime, the system is free of any tachyonic instabilities. The odd winding sectors, hi = 1, 
are exponentially suppressed, and the only significant contributions to the partition function 
occur for hi = 0, ^0 + ^4^4 + ^5^5 = 1 modulo 2 and go -\- 64^4 + 65^5 = 1 modulo 2. 

The pressure receives contributions from the untwisted sector, H — 0, and from the twisted 
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sector H = 1. In the untwisted sector, the result is given by 



R^R^ V ■\ 1 
(2^ 



ai 



E 



(2.21) 



3 • 



"'^ is the number of massless boson/fermion pairs in the untwisted sector of the initially 



supersymmetric A^4 = 4 model. We have 77,^^"*"'=*°^ _ ^^^'2^ where no = 2^ counts the massless 
pairs of the A^4 = 8 model; the factor of 1/2 is due to the orbifolding. 

In the twisted sector we have 



R4R5 1 
(2^ 



(-) 



(e4Q4+e4Q4)94 + (e5Q5+e5Q5)55 



(2.22) 



(2mo + ^o)'i?^ + (2m4 + gA^Rj + {2m, + g,)^R 



3 • 



Tin 



• = 2^/2 is the number of massless boson/fermion pairs in the twisted sector of the 
initially supersymmetric N4 = 4 model. The Qa.Qa.Q^.Q^ appearing in Eq. (2.22) can be 
either zero or identified to the twist charge H = 1. 



Using the definitions of the moduli introduced in Eq. (2.11), the pressure P is taking the 



same form as in Eq. (2.12) with 



E 



(2.23) 



90 + £494 + £595 = 1 mod 2 
90 + £494 + £595 = 1 mod 2 



where the functions ^§09495 given in Eq. (2.15), and the coefficients Ug^g^g^ are similarly 



defined, in terms of supersymmetry breaking R-charges, as in Eq. (2.14), 

2« 



^909495 



^909495 



no ■■ 
no 



2 ' 
1 + 2( 



for Na = 8. 



\(e4Q4+e4Q4)94 + (£5Q5+£5Q5)95 j J\f _ ^ 



(2.24) 



(2.25) 



In the A^4 = 4 cases, the coefficients rig^g^g^ can take negative values as well. The results can 
be generalized to iV4 = 2 models. 
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2.3 The zero temperature limit 



Setting T = 0, or i?o 00, in Eqs (2.10), (2.21) and (2.22), we can obtain the one-loop 



effective potential at zero temperature. It arises from quantum effects due to the spontaneous 
breaking of supersymmetry. For the heterotic models, the effective potential takes the form 



where M is defined in (2.11), and 



v(U) = riioVioiU) + noifoi(W) . 
The coefficients Ug^g^ are determined in terms of the R-charges, 



_A94Q4(s)+35Q5(s) 



(2.26) 



(2.27) 



(2.28) 



and 



^9495 (^) 



7714, rris 



^{2m^ + ~g,)2U-^ + {2m, + ~g,yU 



3 • 



(2.29) 



The type II effective potential takes a form similar to the heterotic one, as in Eq. (2.26), 
where now 



v{U) 



E 



n. 



949: 



95 "^9495 (^) 5 



£494 + esgs = 1 mod 2 
£494 + £595 = 1 mod 2 



with n-g^g^ = rig^g^g^, glveu in (2.25) (or (2.24)). 



(2.30) 



2.4 Non thermal cosmologies 



In the zero temperature limit, the 1-loop effective action takes the form: 



+ 



+ 



V 



(2.31) 



All other moduli can be frozen since they do not appear in the effective potential. More 
precisely, their contributions to the effective potential are exponentially suppressed. Freezing 
further ImS", ImTi and Im?7i, we obtain for the Lagrangian 



(2.32) 
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where 

ReS:=e^^^, T := e^'^^ , W := e^"^" . (2.33) 



It is useful to redefine the fields as follows 



d/V2 -l/V2\ 
I/V2 -I/V2 )\<pr 



(2.34) 



since the field 0_ does not appear in the potential: 



h = -R- -{{d(t)f + {d(l)^f + {d(t)uf) -M^v{U) with ^=^- (2-35) 

We look for homogeneous and isotropic solutions where the metric is of the Friedmann- 
Robert son- Walker form, with vanishing spatial curvature: 

ds^ = -Nitfdt^ + aitfdx'dx' , H = . (2.36) 

Here N is the laps function, a{t) the scale factor and H the Hubble parameter. In the gauge 
choice = 1, the gravitational field equations are 

= ]f + + ^4 + V , (2.37) 

2ii + = -\^^- \^^~ + (2.38) 

Their linear sum is independent of the fields kinetic terms: 

H + 3H^ = V = M^viU) . (2.39) 



Eq. (2.38) follows by differentiating the Friedmann-Hubble Eq. (2.37), once the following 



moduli field equations are satisfied: 

■■ ■ dV 

+ 3if0 = -— = -4M%(W) , (2.40) 
oq) 

dV 

'<Pu + 'iH^u = -w- = -^M^Uv\U), (2.41) 
o<pu 

0_ + 3i/(/)_ = , (2.42) 

where prime derivatives are with respect to U. The last equation can be integrated giving 

1 



2^- a6 



(2.43) 
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where c_ is a positive constant. Eq. (2.41) can be satisfied for a constant W, if there exists 
a solution to 

v'{U) = 0. (2.44) 
We will look for models for which this extremum is a local minimum so that the complex 



structure modulus U is stabilized. The compatibility of Eqs (2.39) and (2.40) requires that 
there exists a constant such that: 

= -4H + ^ . (2.45) 

In [22], it is shown that the solution for = is an attractor. Thus, we concentrate on the 
case C(f, = 0, so that 

M^^ = Mo(^y, (2.46) 
where MoCq is a positive integration constant. 

Higher loop contributions in the effective action are suppressed by powers of Re S" ~ e"'^. In 
particular, they remain suppressed at future times as long as the scale factor a is growing. 
Thus the system remains weakly coupled in the future. If we extrapolate the solution arbi- 
trarily back in time, we encounter in many of these models tachyonic instabilities when the 
supersymmetry breaking scale M becomes of order the string scale, before the Big Bang. 
Eventually perturbation theory breaks down. The analysis of the dynamics of tachyon con- 
densation in this highly stringy regime is beyond the scope of this work. We also point out 
that in certain asymmetric orbifold constructions, despite the breaking of supersymmetry, 
the spectrum is free of tachyons for any value of the radial moduli [12, 19]. It would be 
interesting to analyze cosmological solutions in such type of models. 

The Hubble equation takes the form: 

= -^ + JL where = - c_ > , c = -v{U) x - al^ (2.47) 

a° 5 5 

If c > 0, one has for Cm = 

a{t) = A t^/s where A = 2^/8 Q j . (2.48) 

When the kinetic energy for 0_ is switched on, i.e. when Cm > 0, one has a big bang/big 
crunch cosmology. The solution t{a) is given by 

tia) = ±to [' , 0<a<A where A = ( —] , to = v^c^/^^c"^/^ . 

A/AV1-x10 \CmJ 

(2.49) 
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We investigate whether some of the heterotic and type II models we considered satisfy the 



minimization condition (2.44), which fixes the modulus W, and the positivity of the parameter 



c, Eq. (2.47), which allows for real time solutions: 

Extremum : v'ipl) = , StaUUty : v'iJA) > , Real time : v{U) < . (2.50) 
The shape of the potential as a function of U depends on the R-symmetry breaking charges 



which define the coefficients Ug^g^. Since the functions Vg^g^ in Eq. (2.29) are negative, we 
need some of the Ug^g^ to be positive. 



In the heterotic models, the functions fio and foi defining the effective potential, Eq. (2.27), 
are monotonic functions of U, the former decreasing and the latter increasing, (see Fig. [T]a). 
To have an extremum, nio and noi must be of the same sign. The real time condition requires 



a) ^ 


InU 


b) 


i 


\nU 


/voi 







Figure 1: The functions Vg^g^ versus InU. 

this sign to be positive, which implies that the solution is unstable under fiuctuations of U. 
Stable cosmological solutions exist only at non-zero temperature, as we will see in the next 
section. We may have, however, stationary domain wall solutions. 



We now study the effective potential given in Eq. (2.30) for the different type II cases 



-Case 1 : The potential is proportional to the function fn, which is invariant under 
lA l/U. This implies that a stationary point, which is a minimum, occurs at W = 1, (see 



Fig. [Tp). The conditions (2.50) are simultaneously satisfied if riu > 0. This can be realized 
by choosing = = or H. 

-Case 2 : Since the potential is proportional to the monotonic function foi(Z^), there is 
always a runaway behavior for this model: U i.e. R4 +00 for Uqi > and U +00 
i.e. i?5 +00 for noi < 0. Thus, the dynamics drives the system out of the scope of 
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our analysis, and one should look for a solution in five dimensions, where supersymmetry is 
spontaneously broken by the remaining finite size internal radius. 

-Case 3 : This type II model has an effective potential of the form encountered in the 
heterotic case. 

-Case 1' : The model is supersymmetric with a flat effective potential. 

At finite temperature, we are going to see that the situations in the type II cases also change 
drastically. 

3 Thermal cosmologies 

As we already stated before, we work in a weak coupling regime, with the temperature 
and supersymmetry breaking scales sufficiently below the Hagedorn temperature. It is very 
important to stress here, that the string perturbative computations we use to obtain the 
energy density and pressure are free of any ultraviolet and infrared ambiguities and the 
contributions of moduli that are not participating in the supersymmetry breaking mechanism 
are exponentially suppressed and can be consistently neglected. The solutions we obtain 
describe the cosmological evolution at times much after a Hagedorn phase transition, from 
a high temperature stringy phase. The perturbative analysis is not sufficient to analyze 
the dynamics of the system at the transition, but we expect the solutions to match the 
correct behavior of the system when the Universe has cooled down at temperatures below 
the Hagedorn temperature. In [22], it is explicitly shown that for generic initial conditions 
after the Hagedorn exit, the solution is always attracted to the ones we display below. 



3.1 Equations of motion and thermodynamics 



Using the field redefinitions ( 2.33 )-( 2.34) and the FRW ansatz (2.36), the thermal effective 
action is 



(3.1) 



where the pressure can be expressed in terms of the supersymmetry breaking scale M(0) as 

M 



27r 



u 



(3.2) 
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The gravity equations are obtained by varying with respect to and a, 



iV2 

2 



2m 



+ 



2iV2 



+ 



H + 



HN 



2N^ 
1 



with 



-P-N 



dP 
dN 



P 



1 dP 



(3.3) 



(3.4) 



' j^3-- '2N^^ 2N^^- 2iV2^" ^ 3'^ da 

where in general P can depend on and a. To determine the dependence, we recall in 
which specific frame we computed P: 



N = 27rRoVR^ 



27tRVR^ . 



(3.5) 



where Ri = R2 = R3 = R are the radii of the large three spatial directions (before the large 
volume limit R —>■ +00 is taken). Since P = T^p{u,U) and T is identified with the inverse 
of the laps function, we have 

dP 



dP 



and 



da 



0. 



(3.6) 



It is remarkable that these expressions are identical to the ones derived from thermodynamics. 
We thus show that the variational principle is in perfect agreement with thermodynamics. 



Since P (in the action (3.1 )) and p (in Eq. (3.3)) are scalars under time reparameterizations, 
we can write the gravitational equations in the simple gauge iV = 1 as follows: 

1 



3H' 



2H + 



li2 , li2 



li2 



bl + P: 



1 



P. 



(3.7) 
(3.8) 



Combining Eqs (2.12) and (3.6), we obtain 



p = T^r{u,U) where r = 3]? — uduP ■ 

The equations of motion for the moduli fields can be written as follows: 

■ BP 

(j) + 3H(j)= — =T^{3p{u,U)-r{uM)} , 
dP 

T^V2Udup{uM), 



(3.9) 



bu + 3H<Pu . , 
o(t)u 

L + 3i70_ = , 



(3.10) 

(3.11) 
(3.12) 



where the r.h.s. of Eq. (3.10) follows from Eqs (3.2) and (3.9) 
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The last equation gives (2.43). We would like to find solutions to the remaining system of 



equations, with lA stabilized. Then, Eq. (3.11) amounts to the following algebraic equation 



D{uM) := Udup = 0. 



(3.13) 



requiring that u is also a constant. It follows that the time dependence of p and P arises 
from the pre-factors only. From the relations in (3.2), we have that 



M(0) = uT 



T 
f 



(3.14) 



Instead of solving the scale factor Eq. (3.8), we choose to solve the equation that arises from 



the conservation of the energy-momentum tensor: 



d (\ > 



dt\2 



+ p ) + 3if ^ 



(3.15) 



Using the equations of motion for the scalar fields and (|3.14|), this gives 

-E =^ 
where ao^o is a positive integration constant. 



T 
T 



aT = a(jTo 



(3.16) 



Next, we consider the linear sum of Eqs (3.8) and (3.7) 

1 



H + 3H' 



^{p-p) = \t'{t{uM)-v{uM)} . 



(3.17) 



Using ( 3. 16 ) , the compatibility between this equation and ( 3. 10 ) implies the following thermal 



equation of state: 

p = 5P. (3.18) 

This is one of the main results of this paper. It reminds us of the analogous equation, 
p = 4P, derived in [5,29], when a single modulus was participating in the spontaneous 
supersymmetry breaking mechanism. These results are very suggestive, and we conjecture 
that they will be generalized to the cases when more moduli participate in the supersymmetry 
breaking. When n such fields are involved, we expect the equation of state to take the form: 



p=(3 + n)P. 



(3.19) 



In the two moduli case, which we are considering here, Eq. (3.18 

C{u,U) := {2 + udu)p = 0. 
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) can also be written as 

(3.20) 



As a result, the complex structure ratios {u,U) are determined by the equations D = C = 0. 
It is interesting that along this critical trajectory, complex structure moduli participating in 
the breaking of supersymmetry are stabilized, and thus the cosmology is characterized by a 
single running scale. 



The time dependence of the scale factor is dictated by the Friedmann-Hubble Eq. (3.7). The 
latter takes the form 

3H^ = -\ — ^ where = 6(aoTo)^p(M, W) , = - c_ > . (3.21) 

a° 5 

When Cm = 0, the universe is effectively radiation dominated and a cosmological solution 
exists if the constant p{u,U) is positive: 

a{t) = BVi with B = V2 (y)^^^ • (3.22) 

When Cm is non-trivial, the time t can be expressed as a function of the scale factor as 
follows: 

if > : t{a) = to / , Va > , 

Jo vl + x2 

x^dx 



(3.23) 

ifc^<0: t(a) = ±tol T ^"^ . , 0<a<5, 



a/B Vl - a;2 



where 



B=J^ and to = V3cm\cr\'^^^ . (3.24) 



In the explicit examples presented in the following section, we always find Cr > 0. It would 
be interesting to find if models with "negative effective radiation energy density", Cr/a^, are 
allowed. 

As noticed in [5], the fact that the cosmological evolution we have found behaves effectively 
like a four-dimensional universe filled with thermal radiation is not in contradiction with the 
state equation p = 5P. The reason is that the total energy density and pressure contain the 
"cold" part associated to the kinetic energy of (p. When Cm = 0, one has 0^/2 = p/5, so that 

Ptot = ^(j) + P = 

Ptot = ^0' + P=^p, (3.25) 
in agreement with the expected state equation ptot = 3Ptot- 
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3.2 Some stringy examples 



We examine whether the extremization condition (3.13) and the compatibihty condition 



(3.20) are simultaneously satisfied in the various heterotic and type II models under consid- 



eration: 

Extremum : D{u,U) = , Compatibility : C{u,U) = , Stability : dljp{u,U) < 0. (3.26) 

As in the non-thermal situation, the shape of the potential depends on the R-symmetry 
breaking charges. Their choices determine the coefficients Ug^g^g^ that satisfy — nioo < 

%09495 < "-100- 

However, in the heterotic cases, rim, n-oio and riooi are not totally arbitrary. It is convenient 
to parameterize the a priori allowed models by separating the uiqq states into 4 groups, 
depending on their parity under the operators (— )'^'' and {—)^^, as shown in Table 



3.27 





^loo'Ci states 


^100^2 states 


riioo'Cs states 


^^100(1 - ^1 - 6 - 6) states 




+ 


+ 








+ 




+ 





(3.27) 



We observe that the parameter space of models is the tetrahedron: 

{(^1,6,^3) G [0,1]' such that ^1 + ^2 + ^3 <l| 



which constrains the ratios: 

^010 



^100 



2(^1 + 6)-!, rooi := 



^001 
^100 



2(^1 + ^3)-!, nil := 



(3.28) 



1-2(6 + 6). 

(3.29) 



The conditions (3.26) and the 6 ^ ^3 duality symmetry can be visualized geometrically 
in terms of the tetrahedron representation. Some type II models are also characterized in 
terms of this representation. 

In the large/small u and U regimes, p contains exponentially suppressed contributions that 
we have to neglect by consistency. The dominant contributions take the form of a linear 
sum of a finite number of monomials u°'U^. The (Inw, In W) -plane is divided into 6 sectors 
inside of which a power expansion is defined, (see Fig. [2]). The boundaries of these sectors 
are the lines U = , where uj = —2, 0, 2. We show in the Appendix the power expansion of 
p{u,U) in each sector. More accurate expressions are also given along the lines U = Am^^'° 
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(6) InU (5) 



VI 



(1) 



IV 



(4) 



II 



(2) 



/// 



(3) 



Figure 2: The (liiu,liiU)-complex plane can be divided in 6 sectors, VI, separated by 6 edges of slope 

—2, or 2. The power expansions of p{u,lA), which are listed in the Appendix are well defined in each sector. 

that are useful to connect the power expansions on each sides. D{u,U) = and C{u,U) = 
define curves which asymptote to the fines U = Au^^'", for wfiicfi we fiave determined tfie 
intercept A as a function of tfie ratios rg^g^g^ . Tfie constraints are simultaneously satisfied if 
tfiese curves meet at a point {ucMc)- Tfien, tfie stability condition, dlip{uc,Uc) < 0, and tfie 
sign of tfie radiation density are determined. Tfie C{u,U) = constraint requires to fiave 
at least one negative rg^g^g^. 

• Our analysis sfiows tfie existence of non-trivial tfiermal cosmological solutions in fieterotic 
models witfi (^4 = Q5, and type II Case 3 models witfi (^4 + Q4 = Q5 + Q5, wfiere tfie 
pressure takes tfie general form 



p{u,U) = nioo[pioo + r{poio + Pool) + pi 
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(3.30) 



witfi — 1 < r < 1. Tfiese models lie along tfie edge ^2 = ^3 = of tfie tetrafiedron. Tfie 
duality symmetry U implies tfiat along tfie axis W = 1, D{u,U) =0. As r varies, we 

find 4 distinct patterns, specified by 



rc3 



2/572(1) 



-0.215 , r,2 = - 



SI 

^5° 



and r^i 



0, 



(3.31) 



wfiere tfie functions in tfie definition of can be found in Eq. (A. 12). 



- It turns out tfiat tfiere are no cosmological solutions witfi constant u and U wfien r < r^s 
or r > 0, since tfien C{u,U) 7^ every wfiere. 

- A stable cosmological solution exists wfien < r < rc2, witfi Cr > 0, (see Fig. [3]a). It 
corresponds to a global minimum of tfie tfiermal effective potential. 
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Figure 3: For a pressure of the form (3.30), the curves defined by C{u,U) = (straight lines) and 
D{u,lA) = (dashed lines) are represented. When r^z < r < rc2 (Fig. a), there is a stable cosmological 
solution (uc,Uc = !)• It corresponds to a global minimum of the thermal effective potential. When rc2 < r < 
(Fig. b), there is a stable cosmological solution {uc,Uc = 1) that corresponds to a local minimum. Two run 
away behaviors that bring the system to five dimensions are also allowed. 



- When < r < 0, a stable cosmological solution with Uc = 1 still exists, with Cr > 0, 
but this corresponds to a local minimum of the thermal effective potential. Actually, new 
branches of the locus D = are present and converge exponentially towards the curve C = 0, 
(see Fig. Formally, their common asymptotes define flat directions. However, since these 
"solutions" imply (or i?4) to be very large, they are out of the scope of our analysis. They 
are better understood in terms of runaway behaviors that decompactify the system to five 
dimensions, where supersymmetry is broken by the remaining finite size internal radius R4 
(or Rr^) and thermal effectsP 



(3.32) 



For the heterotic model with the choice Q4 = Q5 = 'y + 7', one has 

nioo = riui = 2^ X 504 , 

noio = riooi = 2^ [ [2]x,,, + [6]t6 + [120 - 128]es + [120 - 128]^^ ] = -2^ x 8 , 
so that < r = —1/63 < Vci- In this specific case, the cosmological solution corresponds 



to {UcMc 



(1.649, 1), where ~ 0.0708 x 6(aoTc 



oj 



Another model considered in [5] is based on the heterotic T^/Z2 orbifold, with non-Abelian 
gauge group EqX EjX SU{2). In that case, the parameter r satisfies < r = —1/127 < 
and the corresponding cosmological solution fixes {ucMc) — (1.996, 1), with Cr — 0.0762 x 
6(aoro)l 



A similar runaway behavior can be realized in type II Case 2 models, with generic operator Qs- They 
involve a pressure p{u,U) = ?iioo[Pioo + ''Pooili where — 1 < r < 1. Their representative points in the 
tetrahedron satisfy ^1 = £,3 = ^ — £,2, where < ^2 ^ |- When r varies, there is a phase where the curves 
D = and C = are non-trivial and asymptotic to one another. The situation is similar to what is observed 
in the lower half plane of Fig. 
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• We can treat in a similar way type II Case 1 ' models with arbitrary Qi + Q^. The pressure 
is of the form 

p{uM) = '^loofPioo + rpul] , (3.33) 

where — 1 < r < 1. This class of models belongs to a segment in the interior of the 
tetrahedron, ^2 = ^2 = ^ — ^i? < < |. The symmetry U U^^ implies D{u,U) = 
along the axis U = 1. Also, p{u,U) is constant in sectors /// and IV (and their common 
edge (4)), implying that D{u,U) is vanishing. For r < 0, there is no other solutions to D = 0. 
When r varies, the set of solutions to C{u,U) = is divided into 3 classes characterized by 



-0.77 



c3 



n 



-0.215. 



(3.34) 



- For r < rc4, the right boundary of the locus C < is asymptotic to the edges (3) and (5), 
where D is not vanishing yet, (see Fig. |4]a). The only solution to D = C = arises at W = 1, 
but the corresponding cosmological evolution is unstable under small fluctuations of U. 

- For < r < r^g, a cosmological solution with constant u and lA exists, (see Fig. |4^), but 
is again unstable under small fluctuations of lA. 

- For r^3 < r, one has C{u,U) > everywhere: There is no cosmological solution with 
constant complex structures. 
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Figure 4: For a pressure of the form (3.33), the curves defined by C{u,lA) = (straight lines) are repre- 
sented. The locus D = is composed of the axis 14 = 1, and the sectors III and IV , (see Fig. When 
f < Tc4: (Fig- 0,) or < r < r^g (Fig. h), there is an unstable cosmological solution (uc,hic = 1)- 



• The last class of systems we analyze corresponds to type II Case 1 models with arbitrary 
Qa + Q^- The pressure is 

p{uM) = '^loofPioo + rpoii] , (3.35) 

where — 1 < r < 1. In contrast to the previous cases, these models are not represented by 
points in the tetrahedron. For r < 0, the only solutions to D{u,U) = are along the axis 
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U = 1. The set of solutions to C{u,U) = is characterized by the point r^g of Eq. (3.34) 
and = 0. 

- For r < r^3, the axis W = 1 is entirely inside the region C{u,U) < that has two distinct 
boundaries, (see Fig. [sja), and so there is no cosmological solution with constant u and U. 

- For r^3 < r < 0, the region C{u,U) < has now a connected boundary, (see Fig. [5^). 
The latter crosses the axis W = 1, so that a cosmological solution {ucMc = 1) exists. It is 
however unstable under small fluctuations of U. 

- For < r, one has C{u,U) > everywhere: There is no cosmological solution with constant 
u and U. 
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Figure 5: For a pressure of the form (3.35), the curves defined by C{u,U) = (straight lines) and 
D{u,lA) — (dashed lines) are represented. When r < r^g (Fig. a), there is no cosmological solution 
with constant complex structures. When r'^^ < r < (Fig. h), there is an unstable cosmological solution 
{UcMc = !)• 



3.3 Non-thermal vs thermal cosmologies 



The cosmological evolutions found in the non-thermal and thermal cases, Eqs (2.47) and 



(3.21), have drastically different properties we would like to comment on. 

First, the contributions of the kinetic energy of the field 0_ to the "effective" Friedmann- 
Hubble equation are of opposite sign. This phenomenon is due to the different compatibility 
relations (f) = —AH for the non-thermal cosmology and (f) = —H for the thermal one, once 



they are inserted in the respective Friedmann-Hubble Eqs (2.37) and (3.7). This is also the 
reason why these equations have and monomial contributions, respectively. The 

intermediate cosmological region we are considering is described by the thermal cosmologies. 
At late times however, the temperature, being proportional to the inverse of the scale factor, 
tends to zero, and a paradox seems to arise as the cosmological evolution is never described 
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by the non-thermal solutions. 

To better understand this point, we rewrite the pressure and energy density as follows, (see 



the relations (3.2) and ( 3.9[ )): 



P 



-M* , p = {v,i, + uduV.t,) where v,t,{u, U) 



(3.36) 



When the pressure is of the general form (A.l), one can use the expansion valid for m 3> 1 



{U fixed), Eqs (A.ll) and (A.15), to decompose v^i, into two pieces: 



U) = - {noiokriK-') + nooikriK) + n^iikl" iU)] + v,^{uM) , (3.37) 



where 



X {('^loo + ^iii)5'4 + (nolo + ^001 + ^oiO^'^ + 



(3.38) 



up to exponentially suppressed terms as m — +oo. Clearly, Vti,{u,U) ^ in this limit. Since 



n 



909495 



n 



9495 5 



the M-independent terms of the r.h.s. of Eq. (3.37) are equal to the zero 



temperature effective potential v{U) found for the heterotic or type II cases, Eqs (2.27) and 
(2.30). Thus, P and p are only converging to their T = counterparts when u +oo. 
However, the thermal cosmologies we have considered have stabilized u = Uc, implying 



that the finite temperature corrections 'Oth(u,W), even if the temperature is small, are never 
negligible. In fact, the condition u = M{t)/T{t) +oo for the thermal system to be 
correctly approximated by the non-thermal one implies that the thermal corrections should 
be screened by radiative corrections, and not that they would necessarily vanish. Explicit 
cosmological evolutions with u +oo are analyzed in [22]. They describe Big Crunch 
cosmologies, where the temperature is proportional to 1/a. So, T is large in absolute value 
(but still negligible compared to M). 

Finally, we note that in the present work, each time a stable cosmological solution at finite 
temperature is found, its counterpart at T = is unstable, and vice versa. As explained 
in [22], this is a consequence of the fact that for arbitrary initial conditions, the cosmological 
evolutions are always converging to an attractor. Depending on the model, the latter can 
be the radiation dominated universes studied in [5] and the present paper, or Big Crunch 
cosmologies where the thermal effects are screened by radiative corrections. 
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4 Conclusions 



The main result of this work is to show the existence of the critical stringy cosmologies 
of [5], even in the cases where more than one modulus participate in the supersymmetry 
breaking mechanism. They correspond to a radiation dominated era with constant complex 
structures. A thorough analysis was done for several string models with N4 > 2 initial 
supersymmetry. Depending on the pattern of supersymmetry breaking, the critical values 
of the thermal effective potential for the complex structure moduli are either minima or 
maxima. Run away behaviors that bring the system to higher dimensions can also occur. 
When the radiation era is stable, we explicitly show that even if the temperature tends to 
zero at late times, thermal corrections to the dynamics are never negligible. 

Our approach was to separate the cosmological evolution in distinct regions, according to 
the value of the temperature. The stringy non-geometrical region, where the temperature 
is of order the Hagedorn temperature, is much harder to understand. Stringy phenomena 
occur where conventional field theoretic notions concerning the geometry and topology are 
breaking down. Some interesting proposals to understand this early time region have been 
put forward in [11, 12]. The intermediate region, just after the Hagedorn era is shown 
to be under control. The free energy is free of any infrared and ultraviolet ambiguities, 
allowing us to follow the backreaction on the geometry and determine the time evolution of 
several moduli fields. We show that the only relevant moduli are the ones participating in 
the supersymmetry breaking mechanism. The others are either frozen by receiving a soft- 
breaking mass, or they remain flat directions with exponentially suppressed contribution to 
the free energy. 

A very interesting result is that the critical cosmological solutions are governed by a higher 
dimensional thermal state equation, p = — 1) P, or in the four-dimensional effective de- 
scription, this equation becomes p = {3 + n) P where n is the number of moduli participating 
in the breaking of supersymmetry. In [5,29], we analyzed the case n — 1. Here, we generalize 
the result to the n — 2 case. The overall scale evolves in time such that the state equation 
is valid, while complex structure-like ratios of scales are frozen. 

Although our work covers adequately and unambiguously the intermediate cosmological re- 
gion, it is incomplete to describe the very early Hagedorn era, where non-geometric stringy 
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phenomena are important. It is incomplete as well to describe relatively late time phenom- 
ena, like the radiative breaking of the electroweak gauge symmetry and QCD confinement, 
where non-perturbative transmutation scales, like for instance Aqcd, are relevant. The 
extension to the late time era requires to consider at least models with A^4 = 1 initial super- 
symmetry. Progress concerning the very early era can be made provided that we understand 
better the resolution of the Hagedorn instabilities and the stringy non-geometrical structure 
of the early universe. 
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Appendix A 

We list in this Appendix the large and small complex structure expansions of an arbitrary 



linear combination of functions Pgog495 defined in Eq. (2.15) 



p{u,U) = nioo(pioo + '^oioPoio + ^ooiPooi + riiiPiu + ronPoii) ■ (A.l) 

For fixed u and A, we determine the power expansion along the lines \nU = culnw -|- InA 
when \nu and/or \nU are large, by making an extensive use of the approximations 

>^ 1 _ 37r 1 1 

t^((2m + l)% + af " 16^7^^'"' 

>^ 1 _ 37r 1 1 ^^■'^> 

^ ((2m)2x + a)^ ~ 16 ^ 7^ ^ " ' ' 
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where a > 0, x — > 0+ and the dots stand for /x'^^^ j terms. Similarly, we use 

^ 1 _ 2 1 1 

t^((2m + l)2x + af/^"3^7^^'"' 
^ 1 _ 2 1 1 

^ ((2m)2x + af^ ~ 3 ^ 7^ ^ ■ ■ ■ ' 

where again we neglect exponentially suppressed terms. The (Inw, lnW)-plane can then be 
divided in 6 sectors, /, //, . . ., VI, where the expansions are independent of uj (and A). The 
boundaries of these sectors are the lines whose slopes are uj = —2,0 or 2, (see Fig. [2]). In 
each sector, we find, 



pi 

nioQ 



pii^ 

nioo 



piv^ 



pV^ 

^100 



M-2 

+u'U'/^i{l + row + rooi)St + (nn + rou)Sl) , 

u^U~^ rooi5° 
+u-'U-'/^ {SI + rooiSD 

WU^'^ ((1 + rooi)^| + (roio + nn + roii)^4°) , 

+u^W ((roio + roii)55" + rooi5f) 

+ (1 + riii)Sl + (roio + rooi + rQii)Sl , 

+u^U-^ ((rooi + roii)^5° + ^oio^f) 

+ (1 + riii)^^ + (roio + rooi + roii)^^ , 

u^W roio^e" 

WU-^/^ ((1 + roio)^: + (rooi + rm + roii)54°) 
u-^ ^° 

WW^I^ {Sl + rowSl) 

WU-^'^ ((1 + roio + rooi)^4' + (^m + ^oii)^, 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 
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where we have defined 
1 

(2m + 1)6 240 



^ 7r3 ^ (2m + 1) 

m 



2_ Stt ^ 1 ^_93_C(5) cje^lSTTV — = A ^ 

5 I6^|2m + l|5 128 7r2 ' ^ 7r3 16 ^ |2m|5 128 7r2 ' (^-^O) 

_ 2 37r 2 1 _ TT^ 2 37r 2 ^ 1 _ tt^ 

^ ~ 7r3 16 3 ^ (2m + 1)^ ~ 192 ' ^"^163^ (2m)4 ~ 2880 ' 



The previous sectors are separated by edges, (1), . . . , (6), in the neighborhood of which some 
terms we neglected in the interior of the adjacent sectors are not exponentially suppressed 
anymore. Along these edges, one has U — Xu^'^''^, where A ~ 1, and 

.(1) 



riioo 

(2) 



P 



riioo 

(3) 



P 



(4) 



P 



nioo 

(5) 



+u' iflJ,{X) + roiof^J,{X-') + rooi/572(A) + (rm + roii)f^J,iX)) , X^U, 

= u^U^/^ ((1 + rooi)-5| + (roio + rm + roii)-5|) 
+^-2 (5r3-(A) + rooi^3"(A)) , A = 

= u^U-^rooiS^ 

+^5/2(A) + {row + rou)hl%iX) + r^oiK/^{X) + ri^iK/2W > A = Uu\ 

= (l + riii)S'^ + (roio + rooi+roii)S'^ 

(roio^3"(A-') + rooiA;r(A) + roiiA;3-(A)) , X^U, 



P 



(6) 



+^5/2(A"') + r^ioKniX-') + (rooi + roii)/i^%(A-i) + niihl%{X-^) , A = Wm'^, 



P 



^5/2(A ^) -r ' 010"'5/2V^^ ; "r V 001 "r ' 01i;"'5/2V^^ 7 "r ' lll"'5/2 

((7|"(A-i) +roioC(A-')) , X^Uu\ 



(A.ll) 



where we have introduced the functions 



•^5/21 ; 3 iQ ((2m)2A+(2n)2A-i)5/2' 

(m,n)^(0,0) ' ' \ J ) 
toe f\N _ A \ 

hi2\^) -^^1^ //2m + 1)2A + (2n)2A-i)5/2 ' (^-12) 

foo / ^ A 1^ V ^ 

■^5/21^; 7r3 16 £^ ((2m + 1)2A + (2n + l)2A-i)5/2 ' 
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2 ^ \ 2 ^ 1 

dl^W = ^Y1 ((2m)2A + (2n + 1)2)3 ' ^3 W = ^Yl ((2m + 1)2A + (2n)2)3 ' (^-13) 



m,n 



7r3 16 ^ ((2m + 1)2A + (2n)2)5/2 



m,n 



5/2 V ^ 7r3 16 ^ ((2m)2A + (2n + 1)2)5/2 ' 

2 ^TT 1 (A. 14) 

''5/2l^J 7r3 16 ^ ((2m)2A + (2n)2)5/2 ' 

r = — — a2 V ^ 

^5/2l^; ^3 16 ^ Z_. ((2rri + 1)2A + (2n + 1)2)5/2 ' 

P'=rA^ = — V ^- ^°°rA^ = — V 

7r3 ((2m + 1)2A + (2n)2A-i)3 ' 7r3 ((2m + 1)2A + (2n + l)2A-i)3 

(A.15) 
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